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ABSTRACT 



/K\ms. In the last 30 years, the existence of small and cool magnetic loops (height < 8 Mm, T < 10^^ K) has been proposed and debated 
to explain the increase of the DEM (differential emission measure) towards the chromosphere. 

Methods. We present hydrodynamic simulations of low-lying cool loops to study their conditions of existence and stability, and their 
contribution to the transition region EUV output. 

Results. We find that stable, quasi-static cool loops (with velocities < 1 km/s) can be obtained under different and more realistic 
assumptions on the radiative losses function with respect to previous works. A mixture of the DEMs of these cool loops plus inter- 
mediate loops with temperatures between 10^ and 10* K can reproduce the observed emission of the lower transition region at the 
critical turn-up temperature point (T ~ 2x 10^ K) and below T = 10^ K. 

Key words. Sun: transition region - Sun: UV radiation - Hydrodynamics 



1. Introduction 

In the last 30 years, different theories have been proposed and 
debated to explain the origin of the solar EUV output at tem- 
peratures below 1 MK. The idea that the transition region (here- 
after TR) emission originates from the bases of the hot large- 
scale coronal loops is not confirmed by the measured DEM 
(Differential Emission Measure) which is o rders o f magnitude 
larger than predicted (Gabriel 1976; Athay I ll981h . The excess 
observed emission compared to the predictions of the traditional 
conduction-dominated TR picture has led, in particular, to the 
suggestion that much of the TR plasma is confined in relatively 
small and cool magnetic loops (height < 8 Mm, T < 10^ K), 
strongly connected to the chromosphere, but thermally ins ulated 
from the corona jD owdv et al. 1 119861: iDowdv 1 119931: fpeldman I 
ll983HFeldmairet al. .,2001.1 . 

The existence of this class of magnetic loops, predicted for 
decades, remains far from being established. From the observa- 
tional point of view, they are indeed very difficult to observe. As 
most of their emission is in the UV, traditional chromospheric 
diagnostics from ground, such as H a, only caiTy information 
about their lower boundary. On the other hand, the lack of spec- 
tral information with proper temperature sensitivity and spatial 
resolution from the current space-born EUV instruments capa- 
ble of observing lines forming between 10"* and 10^ K has so far 
prevented any firm conclusions. 

Some recent observations with th e VAULT instrument (V ery 
High Angular Ultraviolet Telescope, Korendvk e et al. IbOOlh in 
the H I Ly-a line at very high spatial resolution, showed loop- 
like structures with estimated temperatures and densities (T - 
10^^-3x10'* K, P = 0.1-0.3 dyne cm-) that could be appropri- 
ate for the low-temperature end of cool loops (Patsourakos et alT] 
l2007tl\burlidas et al. 2010). While such interpretation has been 
debated (lJudge & Centeno ,,2008 ). those observations, together 
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with the persistence of the problem of the excess emission in 
EUV lines below 10^ K, call for further investigations on the 
matter. 

The g eneral properties of s tatic cool loops were first dis- 
cussed bv lHood & Priest I d 19791) an d then studied more in detai l 
from a theoretical point of view by Antiochos & Noci (1986), 
who demonstrated analytically that a mixture of static cool loops 
with different temperatures can account for the observed lower 
transition region DEM. This class of loops has specific char- 
acteristics compared to the better studied coronal loops. While 
large-scale hot coronal loops have dimensions of the order of 
~ 10 Mm, temperatures hi gher than 10^ K, and obey the "static" 
scaling laws described bv lRosner et aTl ( IT9781 hereafter RTV), 
cool loops are low-lying (estimated heights of the order of 1.1- 
5 Mm), nearly isobaric, and with maximum temperature below 
10^ K. Another major difference with respect to classical coronal 
loops is in the role played by the three terms appearing in the hy- 
drodynamical equation of plasma energy conservation: the heat- 
ing rate, the conductive flux, and the radiative losses. Static cool 
loops are in approximate balance between the heating rate and 
radiative losses. Thus, inversely to coronal loops, the conductive 
flux plays a negligible role. 

Callv&Robb (Il99lh . following the study of 
Antiochos & Noci (1986), treated in more detail the sta- 
bility, structure and evolution of these static cool loops, both 
analytically and numerically. They reveal, in particular, that the 
shape of the radiative losses function, A(r), poses restrictive 
conditions on the existence of static cool loops. In the case of 
a realistic radiative losses function with the p resence of the H i 
Ly -g losses peak aroun d T - 2x 10** K (e.g..l Dere et al. Il2009l 
or IColgan et al. I l2008i dotted and long-dashed line in Fig. [1] 
respectively), strictly static cool loops do not exist (see the next 
sections for details). It is important to note that their analysis 
is limited to a minimum temperature of logT = 4.29 K, so 
that they cannot consider the existence of loops with maximum 
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temperature below this value and cannot deal with Ly-a losses 
properly (even in the thin losses case). They obtain static 
cool loops numerically with maximum temperature around 
r = 8 X 10"* K, assuming a dependence of the radiative losses 
function on the temperature below T - 10^ K (dashed line in 
Fig.lTJ. They conclude, however, that this temperature is a factor 
2-3 too low to explain the observed low-TR emission, even 
if they do not show any calculated DEM function to support 
these conclusions. Moreover, they state that "the steep 
radiative losses function even if so beneficial for the existence 
and stability of cool loops r ests on shaky founda tions". 

After the cited work of 'Callv & Robbl ([1991). our paper is 
the first attempt to study numerically the conditions of existence 
of small and cool loops (height < 5 Mm, T < 10^ K) and, when 
they exist, their physical structure and their contribution to the 
TR emission. We show that stable, quasi-static cool loops (with 
plasma velocities < 1 km/s) can be obtained through hydrody- 
namic simulations als o assuming different f orms of A(r), with 
res pect to the \york of ICallv & Robb I (Il99lh . 

iPeter et aTl (|2004 12006') made the first successful attempt to 
reproduce the shape of the DEM curve quantitatively and qual- 
itatively, even at temperatures below log T = 5.3 K. They used 
a forward model in which they synthesize spectra from three- 
dimensional MHD simulations of the whole Sun atmosphere, 
from the chromosphere to the corona. We think that their re- 
sults do not exclude ours. Their simulations may include struc- 
tures that could be related to the kind of loops we are studying. 
However, the loops we describe in this paper would be covered 
by only very few resolution elements in their simulation, and in 
any case resolving the gradients and the dynamics of the rele- 
vant quantities in our loop models would require a much higher 
resolution. Therefore, we regard our study as cornpleme ntary to 
the kind of large-scale simulations bv lPeter et al. I d2006[) . 

In the course of our study of loops with the properties just 
described, we have also found low-lying quasi-static loops with 
temperatures in the rang e 10^ - 10^ K . Following one of the lat- 
est loop classifications (iReale II2OIOI) we should refer to these 
loops also as "cool loops". In order to avoid confusion we will 
refer to them as "intermediat e tempe r ature loops". This class 
of loops was first studied by iFoukal I (Il976h and observed by 
iBrekke et al. I (1199 7). From observations, they are generally de- 
tected in the UV lines and appear to be steady for long times 
even if more variable and dynamic then hot coronal loops, prob- 
ably due to the pre sence of substantial flows (Brekke et al. 1993 
iDi Giorgio et"an i2003; Reale 207oh. Their estimated densities 
are in the range 10^ - 10'" cm"^ (e.g.. iBrown II 1 996h . Even if ob- 
served intermediate temperature loops are often related to flows 
we decided to report however in this paper our findings on their 
existence (at least analytically) in a quasi-static state. Studying 
the role of flows in the conditions of their existence would mean 
introducing a new dimension in the parameter space to be ex- 
plored, and this is beyond the scope of this paper, concentrated 
on quasi-static cool loops. Anyway, the quasi-static intermediate 
temperature loops we find have different physical properties with 
respect to observed loops with that temperature and this may be 
correlated with flows. 

The paper is structured as follows: in Sec. |2] we describe 
the numerical model, and the different radiative losses functions 
adopted in this work are introduced. In Sec. [3] the hydrodynamic 
simulations and the different loops obtained (cool and interme- 
diate temperature loops) are presented and their properties dis- 
cussed and analyzed. Section [331 in particular, is dedicated to 
the calculated DEMs of these loops and to the comparison with 
the observed one. Finally, in the conclusions (Sec.O, the role of 



the cool and intermediate temperature loops in the solar atmo- 
sphere and the possibility to observe them with current instru- 
mentation will is treated. 

2. Numerical modelling and calculations 

The set of hydrodynamic equations for mass, momentum 
and plasma energy conservation for a fully ionized hydro- 
gen plasma have been solved in a unidimensional, magneti- 
cally confined loop of constant cross-section with ARGOS, a 
1-D hydrodyn amic code with the fully adaptive-grid pack age 
PARAMESH (lAntiochos et al. Ill 9991: [MacNeice et al. Il2000h . A 
fully adaptive-grid is necessary to adequately resolve one or 
more evolving regions of steep gradients, for example the thin 
segments of the transition region chromosphere-corona of the 
loop (for coronal loops) or shocks in dynamical simulations. The 
geometry of the loop is determined by the following analytic 
form: 



z(s) - h 



h{\- Vl-(./y)^) 

1- VI - ('L/2r)' ' 



(1) 



where z is the height above the chromospheric footpoints, ,s- is 
the curvilinear coordinate along the field lines, h is the height 
of the loop apex (i.e., Znwx), L is the loop total length, and the 
constant y is defined in terms of h and L as 



L-2h 



J = 



2 VI -4/z/L 



(2) 



This form implies that z{s) = at the footpoints, z{s) has its 
maximum h at the midpoint, s - Q, and dz/ds - +1 when s = 
+L/2. Equations[T]and|2]generate an arched loop of given length 
L and apex height above the chromosphere h. 

The numerical model includes a thick chromosphere at each 
footpoint 26.7 Mm deep, acting as a mass reservoir; the total 
length of the flux tube is therefore L -n 53.4 Mm. The boundary 
conditions imposed are at the two endpoints (rigid wall, fixed 
temperature), which are located at i = ±{L/2 + 26.7) Mm. The 
temperature of this chromosphere is set at T = lO'' K. This is 
something new with respect to previous works, in which the 
chromos pheric temperatu r e was set to T - 2 or 3 x lO'^ K (e.g. 
Antiocho s&Nocil[T986l ICallv & Robbl [19911 ISoadaro et all 
12003 ). Having a chromospheric temperature at 10^ K allows the 
inclusion of the full peak of the H i Ly-a line in the optically thin 
radiative losses function and to explore the temperature range 
(lO'* - 3 X lO'* K) observed by VAULT The inclusion of the Ly-o- 
as an optically thin line is of course only a first step, suitable for 
this exploratory study. Work is in progress to take into account 
the transfer of Ly-a radiation, as well as to include partial plasma 
ionization at the lower temperature end of the simulations. Since 
we take, by definition, the beginning of the chromosphere as the 
level at which the plasma drops below 10"* K, the exact position 
of the top of the chromosphere (s = ±Li/2 at the beginning of the 
simulation) changes during the calculation with the plasma fill- 
ing or evacuating the loop. So, at end of the simulation, we will 
have a new position for the top of the chromosphere s = +Lf/2 
and, consequently, a new value of h - hf, where hf is no longer 
the geometrical parameter defining the shape of the loop, but the 
height of the loop apex above the T = 10"* K level. 

The set of conservation equations for mass, momentum, and 
energy, respectively, in a one-dimensional plasma solved by 
ARGOS is 



d d 
ot as 
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Fig. 1. Radiative losses functions considered in this work. Solid line: power-law segments function, equal to (for logT < 4.95 K) 
and r ' (for log?" > 4.95 K); solid line plus star symbols: a peak mimicking the H Ly-a losses has been added to the previous 
functio n; dashed line: fr om iC ook et al. (1989), but with a dep endence below T — 10^ K; long-dashed line: from iColean et alTI 
(120081) : dotted line: from lOere et al. I (120091) : dot-dashed line: from lDere et al. I (|2009|) without the H contribution. 
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log T > 4.95 K (oc r ), with a and b both positive. In particular: 



where t is the time, p the mass density, v the velocity, f , T and n 
are the gas pressure, temperature, and electron number density, 
respectively, linked by the equation of state P - 2nkT, with k 
the Boltzmann constant; U - 3P/2 is the internal energy, s the 
curvilinear coordinate along the loop, E(s, t) the assumed form 
for the coronal heating rate, n^A(T) the plasma optically thin 
radiative losses, with A(r) the radiative losses function, g\\is) 
the component of the solar gravity along the loop axis, and Fc 
the thermal conductive flux, in CGS units. 



2.1. Radiative losses functions 

The hydrodynamic simulations are performed assuming dififer- 
ent shapes for the radiative losses function, that are shown in 
Fig. [1] Since radiative losses were hard-coded in the original 
version of ARGOS, we modified the code to allow for arbi- 
trary, tabulated A{T). We have verified with some benchmark 
tests, that the modification does not introduce significant over- 
heads in the calculations. First, we use the s ame approximation 
of the radiative losses function adopted by lAntiochos & Noci"! 
(Il986l solid line in Fig. [D Aan in Table [D, who fitted the 
data of lOaetz & Salpeter (1983), assuming that the A(r) can 
be approximated by an increasing power-law of temperature for 
log T < 4.95 K (oc T"), and by a decreasing power-law for 



Aam(T) = 



7.94 X 10-^^r^ for logr < 4.95 K 
5.62 X lO '^r ' for \ogT > 4.95 K 



(7) 



This approximation does not consider the peak due to the H Ly- 
a radiative losses at around T = 2 x 10^ K. In order to study its 
effect on the simulations, we added to the function Aan a peak 
mimicking the H Ly-a losses (solid line plus cross symbols in 
FigH] Aanlw in TablelTJ. The peak has approximately the same 
shape of the one in the realistic radiative losses function Acea, 
that we are going to describe. It was built by taking the two lines 
parallel to the ascending and to the descending branches of the 
Ly-ff peak present in the function Ac™, and making them pass 
through the point with abscissa log T - 4.21 K, corresponding 
to the peak top value. 

The long-dashed line and the dotted line in Fig. [1] are two 
re alistic radiativ e losses functions. The first is from the work 
of lCoIpan et"an (|2008. Acea in Table [1) and the other was ob- 
tained using the "CHIANTI" database, version 6 (iDere et al. I 
20091 An«„ in Table [B, adopting the ionization equilibrium of 
Brvans et al. I (l2009l) . and a coronal mixture of elements. We 
choose these two functions because they are the most recent and 
are based on the latest atomic data available. A discussion on the 
differen ces, quite sig nificant, between the two functions can be 
found in lDere et al.T(l2009i) . 

The last radiative losses function used in this paper has been 
obtained by removing the contribution of hydrogen from the 
function Apea (dot-dashed line Fig.lT] Aoea-H in Table [TJ. Even 
if the hydrogen losses have been subtracted, a smaller peak re- 
mains at around T - 2 x 10'* K due to other elements losses 
(mainly helium). 

In the next section, the results from simulations using these 
different radiative losses functions will be introduced. The list of 
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simulated loops, together with the relevant parameters, is given 
in Table E 

3. Results and discussion 

To derive the i nitial equilibrium, we followed the approach of 
ISpadaro et al. 1 (l2003l) . i.e., we simply started with a discontin- 
uous temperature profile in which T - Ti (with T; a constant 
value) in the corona and T - 10"^ K in the chromosphere. The 
density was constant in the corona, n = rii, and increased expo- 
nentially with depth in the chromosphere at the appropriate scale 
height. Hence, the plasma was in approximate force balance to 
start with, but not in thermal equilibrium. During the simula- 
tion, the loop evolves towards a stable quasi-static solution of 
the hydrodynamic equations (Eqs.[3]-|5]l, a solution, that is, that 
approximately fulfills the static and stationary versions of those 
equations: 

d 

—P = -mungnis), (8) 
as 

|-Fc = Ei,-n^A(T), (9) 
as 

where we used p = mun for a fully ionized hydrogen plasma 
with ;«H the hydrogen mass. We consider a loop in a quasi-static 
equilibrium state when the plasma velocities are lower than 1- 
2 km/s. With this definition, we find that in all our simulations 
the thermodynamic parameters of a loop in quasi-static equilib- 
rium oscillate around the equilibrium value with an amplitude 
much smaller than 1%. We took this quasi-static equilibrium 
state as a result or, in some cases (see below), as an initial equi- 
librium state (f = s) for a new simulation with different param- 
eters (e.g.: different heating rate). 

In all the simulations we assume spatially uniform and tem- 
porally constant heating rate per unit volume, E(s, t) - E\^, and 
we ensure that they last long enough for the end time to be much 
larger than the loop radiative cooling time. 

3.1. Description of results 

We ran numerous simulations, extensively exploring the parame- 
ter space, under different initial conditions. We list in Tab.[T]and 
discuss below a representative selection of the loops in quasi- 
static equilibrium we have obtained. Loops 3-4, 7-13, and 17- 
23 are obtained by starting the simulation with a loop in ap- 
proximate force balance but not in thermal equilibrium as de- 
scribed previously, using different combinations of T, and tii that 
we thought be appropriate for cool loops. All these simulations 
have an initial half length, L-Jl - 6.5 Mm and an initial height, 
hi = 0.62 Mm. Once we obtained two quasi-static cool loops 
(loops 3 and 4), we used them as initial equilibrium states to 
perform new simulations. In this way, starting in particular from 
loop 3 and changing some parameters we have obtained all the 
other cool loops 1-2, 5-6, 14-16, and 24-26. Loops 1 and 2 have 
been obtained by changing the heating rate and the initial height, 
hi - 0.7 Mm. Loops 5 and 6 have been obtained by chang- 
ing again the heating rate and the initial height but to the value 
hi - 0.54 Mm. Finally, loops 14-16 and 24-26 have again differ- 
ent heating rates but also different radiative losses functions. 

3.1 .1 . Loops from idealized A(r) 

The first group of simulations we have done assumes as a ra- 
diative losses function Aan, in order to test the analytical re- 



sults of lAntiochos &Nocil(II986h . We are able to obtain quasi- 
static cool loops with maximum temperature between ~ 2.6 and 
8.3 X lO'* K, using E]^ in the range 4 - 7.5 x lO"'* ergs cm"-' s"'. 
In Table[T]we report only some example of the cool loops found 
(loops 1- 6). These loops hav e the properties analytically pre- 
dicted by lAntiochos & Nocil ([l986): they are small (L/2 = 
0.7 - 2.5 Mm and h = 0.005 - 0.04 Mm), nearly isobaric, and 
in approximate balance between the heating rate and radiative 
losses, Eh = n^A{T) » |VFc|. Also the low pressure values, 
compare d to coronal hot loops, were analytically predicted by 
Antiocho s & Noci 1 ([T986h . According to their calculations, the 
pressure of a cool loop has to be at least one and one-half or- 
der of magnitude smaller then that of a hot loop with maximum 
temperature bigger than 10^ K. This result is confirmed by our 
calculations. 

In Figs. 12] and [3] we plot the behavior of the loop parameters 
as well as of the terms of the energy equation |9] for five loops 
chosen as examples for each group of loops obtained with a dif- 
ferent A(T) (loops 4, 10, and 15 from top to bottom in Fig. |2] 
and loops 17 and 25, from top to bottom in Fig. [3]). The left pan- 
els show the temperature (solid line) and the pressure (dashed 
line) profiles as a function of the curvilinear coordinate, s, while 
the right panels show the radiative losses energy term, n^A(T) 
(crosses), the heating rate, E]^ (solid line), and the divergence of 
the conductive flux, VFc (asterisks), as a function of the temper- 
ature. The plots show the variation of loop temperature, density, 
and energy terms along the loop at a particular instant (corre- 
sponding to the end time of the simulation). Since the loops are 
still in a quasi-static equilibrium state, with the plasma subject 
to motions (oscillations, for instance) of small, but not negligi- 
ble amplitude, those profiles do not appear exactly symmetric 
around the loop apex {s = 0). Averaging profiles over an appro- 
priate number of time steps would cancel out such small varia- 
tions, thus yielding nearly symmetric profiles. 

The maximum temperature of loop 4 (top left panel of 
Fig. |2]l, chosen as the example for the first group of loops 1 -6, is 
8.3 X 10"* K and the pressure is constant along the loop (within 
1% above the chromosphere). The terms n^A{T) and E\^ are in 
approximate balance, while the divergence of the conductive flux 
is only a small term (right panel). 

Using the same radiative losses function, we obtain also 
quasi-static intermediate temperature loops with maximum tem- 
perature r ~ 8 X 10^ K (loops 7-13 in Table [B. For loop 10, we 
show in the middle panels of Fig.|2]the behavior of the tempera- 
ture and the pressure as a function of s (left) and of the terms of 
the energy equation as a function of the temperature (right). The 
two terms competing in the energy balance are the divergence 
of the conductive flux (asterisks) and the radiative losses term 
(crosses) as in hot coronal loops. 

We explore the possibility of obtaining stable cool loops also 
with the function Aanlvu- Using the cool loop 3 of Table [T] as 
a starting solution, we made a simulation by changing the ra- 
diative losses function from A^^v to AANLya- The temperature of 
loop 3 drops below T = 2 x 10"* K (the position of the Ly-a 
peak) and a quasi-static, stable, cool loop of maximum temper- 
ature 1.2 X 10"* K forms (loop 14 in Table [1]). The temperature 
of this loop is comparable with the loop-like structures observed 
by VAULT but the pressure is ten times lower. Loop 14 is also 
lower and shorter with respect to loop 3. In order to obtain hotter 
loops (but still cool, with T < 10^ K) from loop 3 by chang- 
ing only E]^ and using Aanlw, we have to raise its value from 
5.5 X 10""* ergs cm"-' s"' to 5.8 x 10""* ergs cm"-' s"' . The increase 
of E]^ should compensate the higher radiative losses (because of 
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Fig. 2. Top panel, left: temperature (solid line) and pressure (dashed line) as a function of the curvilinear coordinate along the field 
lines, s. Right: divergence of the conductive flux (asterisks), radiative losses (crosses) and constant heating rate, E]^ (solid line), as 
a function of the temperature, for loop 4. Middle and bottom panels: as in the top panels for loop 10 and 15, respectively. 



the introduction of the Ly-a peak) and allow the loop to keep the 
initial equilibrium without dropping to low temperatures. We ob- 
tain, indeed, a quasi-static cool loop with maximum temperature 
r = 7.4 X lO-* K (loop 15 in Table[B. Loop 16 of Table[T]is in- 
stead obtained with = 5.9 x lO"'* ergs cm"-* s"' and Aanlvo- 
Using - 6 X 10 ergs cm"^ s the loop reaches a maxi- 
mum temperature much higher than 10^ K. So, loop 15 and 16 
are the two loops with temperatures higher than the temperature 
of the Ly-a peak that can be obtained from loop 3 changing only 
E^i and the radiative losses function to Aanlvu- We are then able 
to "overcome" the Ly-a peak, whos e presence was thoug ht to 
prevent the existence of cool loops jCallv &Robbl[T99l . The 
three cool loops 14-16 have the same characteristics already de- 
scribed for loops 1-6. In the bottom panels of Fig. |2l we show 
the behavior of the temperature and the pressure as a function 
of s (left) and of the terms of the energy equation as a function 
of the temperature (right) for loop 15. The temperature of the 

loop starts to increase at around s = -1.1 Mm up to i 1 Mm, 

reaching the value of ~ 1 .2 x 10"* K and then increases rapidly till 
the maximum value of ~ 7.4 x lO'* K. Loop 16 has the same be- 
havior. The values of L/2 in Table [1] for these two loops include 
the piece where the temperature rises slowly. The terms compet- 



ing in the energy balance are the radiative losses term (crosses) 
and the background heating (solid line) as for loops 1-6. 

The intermediate temperature loops we have found are 
longer than cool loops. Since the loop height is linked to the loop 
length, lower loops are also shorter. In the next sections, we are 
going to show that quasi-static cool loops need to be low-lying 
to satisfy the quasi-static equations with a negligible conductive 
flux. Intermediate temperature loops, instead, do not need to be 
as short as cool loops since their length is linked to the pressure 
and to the heating rate by RTV-like scaling laws. 

3.1.2. Loops from realistic, optically thin A{T) 

Using the function of radiative losses Acea, we obtain only 
quasi-static loops with peak temperature higher than 10^ K 
(loops 18-23 in Table [TJ. Loop 17 has been obtained with the 
realistic radiative losses function Aoea but has the same proper- 
ties of loop 18. For loop 17, we show in Fig. [3]the behavior of 
the temperature and the pressure as a function of s (left panel) 
and of the terms of the energy equation as a function of the tem- 
perature (right panel). The energy balance terms behave as for 
intermediate temperature loops 7-13. 
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Fig. 3. As in Fig.|2]for loop 17 (top panels) and loop 25 (bottom panels). 
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Using Aoea-H and following the same method used in the 
case of Aanlxu, we obtain the cool loops 24-26 of Table [1] The 
maximum temperature of loop 24 is lower than the temperature 
at which Aoea-H shows the small peak (at T ~ 2 x lO'* K) that 
remains even if the hydrogen losses have been subtracted. Loops 
25 and 26 have instead a higher temperature. For loop 25, in the 
bottom panels of Fig.[3]we show the behavior of the temperature 
and the pressure as a function of s (left panel) and of the terms 



of the energy equation as a function of the temperature (right 
panel). As for loops 15 and 16, the temperature of loop 25 starts 
to increase at around s = -5.3 Mm up to s 0.8 Mm, reach- 
ing the value of ~ 1.2 x lO'* K and then increases quickly till the 
maximum value of ~ 4 x 10"^ K. Loop 26 follows the same be- 
havior These two loops are longer with respect to loops 15 and 
16 if we consider the piece where the temperature rises slowly. 
We think that this difference is due to the different slopes of the 



Table 1. Loops parameters. 
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Fig. 5. Evolution of the maximum temperature of loop 11 in Tab. [T] during the hydrodynamic simulation. The horizontal solid line 
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A(r) used in the two cases for T < 1 .2x10'' K. So, considering a 
total height of 0.29 Mm for loop 25, only the top part with length 
of ~ 0.005 Mm is hot. This loop is really at the limit because of 
its dimensions and its shape may depend from the details of the 
boundary conditions and the radiative transfer as we suggested 
previously (the shape of A(T) around T - lO'* K). In order to 
observe these kind of loops we would need at least two different 
spectral lines formed at temperatures between lO'* and 2 x 10"* K 
in order to resolve the loop in its te mperature extension and a 
very high resolution. As predicted bv Antiochos & Noci I (Il986h . 
also for loops 24-26, the terms competing in the energy balance 
are only the radiative losses term (crosses) and the background 
heating (solid line). 



3.1.3. Relations between loop parameters, and scaling laws 

In Fig.|4]we show the relations between the thermodynamic pa- 
rameters (P, Tmax and L/2) and the heating rate for the loops 
in Table [T] (loops 1-6 are represented by diamonds, loops 7-13 
by crosses, loops 14-16 by triangles, loops 17-23 by asterisks, 
and loops 24-26 by squares). The solid lines in the lower panels 
of Fig. |4] represent the RTV scaling laws for coronal loops for 
different values of L/2, while the dotted line (lower-right panel) 
represents the law - P^. The pressure of all the cool loops 
with T < 10^ K is proportional to the square root of and 
it is independent from their length and maximum temperature. 
Intermediate temperature loops 7-13 and 17-23 (with different 
proportional coefficients) obey the RTV scaling laws for coronal 
loops. This is a somewhat surprising result since their tempera- 
tures are lower than 10* K and observed interm ediate tempera - 
ture loops do not o bey the corona l scalin g laws dBrown 
The model used by iRosner et al. I (Il978h to derive the relation- 
ships between coronal temperature, pressure, length and heating 
in coronal loops assumes a global energy balance between the 
heating and the radiative losses in static conditions. The majority 
of observed intermediate temperature loops reveal temperatures 
and densities that lie outside the static model tracks and, as dif- 
ferent authors have shown (e.g., Dere 1982; Durrant & Brown 
11989 ). the static model does not seem to accurately predict their 
physical conditions. The difference between our loops and the 
observed ones is in the pressure. We obtain intermediate tem- 
perature loops with pressures that are 1-2 orders of magnitudes 
lower t han mea sured in observed loops with the same temper- 
atures (lBrownlll996h . Loops 7-13 and 17-23 actually have the 
right pressures to fall on the scaling laws lines. 



3.2. On the conditions of existence of cool loops 

The quasi-st atic, stable cool loops w ith the characteristics 
predicted by l Antiochos & Nocil (11986 1 found are un e xpecte d 
according to the considerations of Callv & Robb I (|1991|) . 
Assuming a form of oc (n/noY, with «o the density at the 
loop base, the conditions of existence and stability of cool loops 
derived in the mentioned paper are not valid in the case of v = 
(constant heating per unit volume) combined with a - 2. This 
combination should not allow strictly static solutions. So, it is 
rather remarkable that we obtain quasi-static cool loops with the 
radiative losses functions Aan, A 



AN, J^ANLy a- 



and Aoea-H- The first 
two functions are, indeed, built with a dependence below 
10^ K (excluding the peak for Aanlvo) and, from Fig. [T] it is 
clear that also Aoea-H follows the same behavior (excluding the 
small peak). 



a = 2 is also empiricall y a better approximat ion to realis- 
tic A(r) than a = 3 used bv ICallv & Robb I (Il99lh . They show 
th at if the hydrogen losses are igno red, following the suggestion 
of iMcClvmont & Canfield I d 19831) in order to obtain a better ap- 
proximation to the true radiative losses function, then cool loops 
can exist. They assume, however, that when ignoring the hydro- 
gen losses the peak in the radiative losses function is replaced 
by a steep slope for T < 10^ K. Actually, ignoring the con- 
tribution of the hydrogen in the function Apea, for example, we 
obtain a dependence closer to T^. Looking also at the functions 
Aoea e Acea it IS possiblc to say that the approximation below 
10^ K is closer to the shapes of realistic radi ative losses functions 
than r^, and it is also of ten used (see, e.g.: iRosner et al. 1 fT978t 
lAntiochos &Nocilll986l) . 

More surprising is the existence of cool loops even if the ra- 
diative lo sses funct i on inc ludes a peak around T = 2 x lO'' K. 
For Call v & Robb I ( 1199 lb . its presence prevents the existence 
of cool loops at the conditions we are considering (v - 0) and 
not only. Their hydrodynamical simulations start always from 
an initial uniform cool state (logT = 4.3 K), corresponding to 
the equilibrium solution in the absence of gravity and, they say, 
the presence of the Ly-a plateau (they do not include the whole 
peak) prevents cool solutions, even unstable ones, with the loops 
reaching directly the hot state. All the simulations performed 
with a radiative losses function including the Ly-a peak (AANLm 
and Aoea-H), start by setting the quasi-static cool loop 3 as ini- 
tial equilibrium state (as descri bed in S ec. [3Tb so that we did not 
have the problems faced bv lCallv & Robbl(ll99ll) . 

Moreover, the cool solutions we have found, not only ex- 
ist, but are also stable ( contrarily to previous predictions, e.g., 
lJudge & Centenoll2008h . in the sense that the loop thermody- 
namic properties oscillate around an average value. These os- 
cillations are of small amplitude and do not grow with time. 
Occasionally, however, these oscillations bring the loop to a dif- 
ferent stable state, determined by the particular shape of the 
A(r) adopted. The stability of these loops is therefore related 
to the chosen shape of A(r) and, in particular, to the tempera- 
ture value at which it reaches its maximum or changes slope. In 
the case of Aan, if we change the maximum temperature value to 
a higher one (and increase the value of E]^, we can still obtain 
stable cool loops with a higher maximum temperature. As an 
example, we have included in Table [T] loop 27 obtained by rais- 
ing the temperature value at which Aan peaks from T - \{f '^^ to 
10'' K(Ayi^io5 inTable[TJ, and using L^ = 6x10""* ergs cm"-' s"'. 
The loop, during the simulation, reaches a stable and quasi-static 
configuration in less than one hour with the maximum temper- 
ature oscillating around the value ~ 8.7 x lO"* K. If we change 
^ANHf using the original Aan function, the loop reaches a new 
equilibrium and its maximum temperature increases to a value 
of ~ 8 X 10^ K (loop 11 in Tab. [TJ. This happens because the 
random oscillations of the maximum temperature of the loop 
around its equilibrium value have sufficiently large amplitudes 
to overcome, at some point, the position of the peak of the A(r) 
function, which for Aan is at \ogT break - 4.95 K; if that happens, 
the entire loop suddenly switches to a different equilibrium state, 
corresponding to an intermediate or coronal-type loop. In Fig.|5] 
we show the evolution of the maximum temperature of loop 1 1 . 
In the case of Aoea-H, the maximum temperature that can be 
reached by stable cool loops is related to a change in the slope 
(fl < 2) at logT ~ 4.8 K. If we increase the value of E^ using, 
for example, loop 26 as a starting solution for a new simulation, 
the random oscillations of the maximum temperature around 
the equilibrium value will make it to overcome log?" ~ 4.8 K 
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and the loop will find a hot equilibrium state (ISerio et al. II 198 it 
lLitwin&Rosnerlll993l) . 

The characteristics of the cool loops and, in particular, the 
behavior of the three terms of the energy equation can be 
derived from Eqs. p and |9| In the equations considered by 
lAntiochos & Nocil (Il986h . the conductive flux is zero along the 
loops, so that the energy conservation equation (Eq. |9]l becomes 
Zih = rP'MT). If fl = 2, the radiative losses rp-K(T) are propor- 
tional to (making use of the equation of state) and indepen- 
dent on the temperature, thus the energy balance implies that the 
pressure is constant. This requirement is obviously incompatible 
with pressure variations imposed by the hydrostatic equilibrium 
(Eq.|8]l. So, we need to relax the conditions of zero conductive 
flux along the loop, considering |VFc| <s ish- In this case, we can 
still have some physical solutions. Derivating Eq.|9]with respect 
to s and, using the equation of state, we obtain 



A(T) d 



(10) 



where we have made use of the fact that A{T)l{2kT)^ is con- 
stant, if a = 2. In hydrostatic equilibrium we can easily estimate 
how small the required pressure gradient should be to keep the 
conductive flux negligible with respect to the other terms in the 
energy equation; using Eq.|8]in Eq.[TO] we have: 

ws'^ - ^^^^2^^) wry ^''^ 

where H{T) is the pressure scale height as function of temper- 
ature, defined as H{T) = 2kT jim^g), and the function ri{s) = 
g\\(s)lg{s) is bound in the interval [0,1]. By imposing the ap- 
proximate equality between Zsj, and radiative losses: 

^^-c - £h^. (12) 

An order of magnitude estimate of the loop dimension is there- 
fore: 



d^s h\ds I H{{T)y 



(13) 



where the order of magnitude of the quantities is given by their 
averages along the loop; factors of the order of unity have been 
dropped. Thus, the requirement that |VFc| <s; implies that 
h <K H(T), i.e. the height of the loop must be much smaller 
than the scale height relative to the "average" temperature of the 
loop. Considering, for example, that the pressure scale height at 
r = 6 X lO"* K is 3.6 Mm and IVFcl/^h is ~ 0.03 for our loops, 
then their heights are compatible with this order of magnitude 
estimation. 

In other words, the limitation set to the solutions of Eqs. [8] 
and|9]by the case A(r) oc T^, can be circumvented for quasi- 
static solution, only for low-lying (elongated) loops. Indeed, our 
solutions with temperatures below 10^ K are all a fraction of Mm 
high (the maximum is ~ 0.3 Mm, for loops 24-26 in Tab.[T]i. 

3.3. Calculated DEMs for cool and intermediate temperature 
loops 

The theoretical DEMs for the quasi-sta t ic loo ps we have found 
are computed according to iPeter et al. I (l2006l) . with a tempera- 
ture bin of 0.05 dex on a log T scale, and considering a filling 
factor of 100%: 



T dh 

DEM = — . 

dT 



(14) 



The upper panel of Fig. |6] shows the calculated DEMs ver- 
sus temperature of the quasi-static cool loops 1-6 (solid black 
lines), 14-16 (red), 24-26 (green), and the quasi-static inter- 
mediate temperature loops 17-23 (magenta) of Table [T] In this 
figure and in the next ones, the calculated DEMs are com- 
pared with the observed DEMs of quiet Sun and active re- 
gion (dashed and dotted lines, respectively), derived using the 
Vernazza & Reeves ( 1978) average quiet Sun and active region 
intensities, and produced as part of the Arce tri/Cambridge /NRL 
"CHIANTI" atomic data base collaboration dPere et al. I l2009). 

In the lower panel of Fig. |6] we plot the total theoretical 
DEMs, for each group of loops obtained with a different A(T) 
(distinguished by the different colors). Assuming that the loops 
are equiprobable (uniformly distributed in logr) and with the 
same cross-section, we divide the temperature range in bins of 
amplitude 0.2 dex on a log T scale, and consider for each bin 
a representative loop, i.e. a loop whose maximum temperature 
belongs to that bin (our loops are almost isothermal). The total 
DEMs are obtained by summing the DEMs of these representa- 
tive loops. When more loops have their maximum temperature 
falling in the same bin, we average their DEMs. From Fig.|6] we 
see that the largest contribution to the total DEM for the cool 
loops is given essentially by the peaks at maximum temperature 
and this contribution is determined not by the form of the single 
loop DEMs, but by the distribution of the emission measures 
from loop to loop as predicted by Antiochos & Noci (198^. 
This is a diff'erent behavior with respect to the hot and the in- 
termediate temperature loops. The DEM of a group of coronal 
loops is dominated at all temperatures by the hottest loops, hav- 
ing the largest emission measure. Hence, the DEM obtained by 
summing all the DEMs will closely resemble the DEM of the 
hottest loop. For a more accurate discussion on the shape of the 
resulting DEM from a group of cool and hot loops we refer to the 
calculations of Antiochos & Noci ( 1986). They made assump- 
tions on different variables (for example, the magnetic field) and 
loop distributions that we do not include in our work. 

Using Aan, we can easily obtain loops with maximum tem- 
peratures covering the whole temperature range, from log?" ~ 
4.1 K up to the position of its peak (logT - 4.95 K). Adding 
more DEMs of loops with different maximum temperatures, 
the resulting DEM (black solid line) can be brought to fol- 
low the shape of the observed ones. However, using Aanlw or 
Aoea-H the presence of the Ly-or peak (or a peak in general) at 
\ogT ~ 4.2 K produces a relative minimum in all the DEMs, that 
remains in the total DEM (lower panel of F'ig. |6l green or red 
lines). We note that iMacpherson & Jordan I (|1999) derive quiet 
Sun DEMs that exhibit a shape around logT" ~ 4.2 K that could 
recall the minimum that we find in our results. 

The DEMs of the intermediate temperature loops 20-23 
show a minimum at a different temperature with respect to the 
observed DEMs. We are unable to reproduce the whole observed 
DEMs considering a unique loop of this kind, in analogy with 
the relationship between the DEM of a hotter loop and the ob- 
served DEM of the corona above 10^ K. From Fig. |6] (lower 
panel), it seems more likely that a combination of the cool and 
intermediate temperature loops would assume the right shape to 
reproduce the observed emission of the lower transition region 
at the critical turn-up temperature point (T ~ 2 x 10^ K) and 
below T = lO^* K. In Fig. |7] we show indeed the total DEM 
(solid line) resulting from the combination of the DEMs of the 
cool loops 24-26 plus the intermediate temperature loops 17-23. 
We choose to sum the DEMs of these particular cool and inter- 
mediate temperature loops because they have been obtained by 
using realistic radiative losses functions. Even though the A{T) 
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Fig. 6. Top: calculated DEMs for the quasi-static cool loops 1-6 (solid black lines), 14-16 (red), 24-26 (green) and the intermediate 
temperature loops 17-23 (ma genta) of Tab. [1] c ompared to the DEMs of quiet sun (dashed) and active region (dotted) from the 
"CHIANTI" atomic data base (IDere et al. II2009I) . Bottom: total DEMs for each group of loops shown in the top panel. 



adopted are different, they can still be combined together to ob- 
tain a single DEM. We are legitimated to so because, while for 
the existence of the cool loops the shape of A(r) under 10^ K 
is very important, we have proved instead that the exact form of 
the chromospheric radiative losses has little effect on the coro- 
nal properties of the loops with temperature higher than 10^ K. 
For loop 17, for example, changing the radiative losses function 
from Aoea to hoea-H docs not bring any change in the loop sta- 
bility and thermodynamic parameters. This result extends what it 



is already known for h ot corona l loops (T > 10^ K) and cooling 
coronal loops (e.g.. Reale et al. I [l988: Brown 1993). 

There is a minimum in the total DEM around log?" = 4.8 K 
that is due to the lack of cool loops with that maximum tem- 
perature and Aoea-H- This minimum almost corresponds to the 
maximum of the function hoea-H or better to the point where 
its slope starts to change and we have a < 2 (we have akeady 
considered the consequences in Sec. 13.2b . So, the lack of cool 
loops with maximum temperature around log T - 4.8 K it is 
not due to an incomplete exploration of the parameter space but 
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Fig. 7. Total DEM resulting from the combination of the DEMs of the loops 17-23 and 24-26 (solid line ), compared to the DEMs 
of quiet sun (dashed) and active region (dotted) from the "CHIANTI" atomic data base (iDere et al. II2009I) . 



to the negative slope of Aoea-H that prevents their formation. 
We could use the DEMs of the other cool loops of Table [1] with 
maximum temperature around 6.5 x 10"^ K or higher but they 
have been obtained from idealized A(r) and would not be suit- 
able for a good comparison with the DEMs of the intermediate 
temperature loops 17-23. However, the shape of the averaged 
DEM of the loops 17-23, with a flat minimum and a tail extended 
towards low temperatures, helps filling this gap, improving the 
agreement with the observed DEM. Since we considered a fill- 
ing factor of 100% the total DEM has its highest value. With a 
lower fiUing factor the height of the DEM would be lower. 



4. Conclusions 

We have studied the conditions of existence and stability of cool 
loops with T < 10^ K through hydrodynamic simulations, find- 
ing that it is possible to obtain quasi-static (vel ocities lower tha n 
1 km/s) cool loops, as predicted by Antiochos & Nocil (Il986h . 
stable over hours or more. These loops have been obtained by us- 
ing different dependences of the radiativ e losses function on th e 
temperature with respect to the work of ICallv &Robbl(ll99lh . 
We obtained stable quasi-static loops even in conditions judged 
prohibitive by the same authors on the basis of an analysis of 
strictly static and stationary loop equations. We examine and dis- 
cuss the quasi-static solutions we have found, and show that their 
existence is due indeed to the small departures from static con- 
ditions, i.e. to the presence of a small but non-zero conductive 
flux and velocities, as well as to rather stringent constraints on 
the pressure gradients. In fact, for low-temperature loops, the re- 
quirement of nearly constant pressure implies that these loops 
can exist only if they are limited to small heights above the chro- 
mosphere (a fraction of Mm). We also show that the presence 
of the peak due to the Ly-a losses in the radiative losses func- 
tion does not preclude the existence of cool loops. Moreover, we 
analyze the contributions of cool loops to the TR DEM, show- 



ing that the emission of these kind of loops can account for the 
observed DEM at T < 10^ K, if they were uniformly distributed. 

The cool loops found cannot be related with any of the ob- 
servations present in the literature because of their dimensions 
and especially of their small heights. We find cool loops with 
lengths of 5-10 Mm, but with very low heights (in the range 
10"^ - 10"' Mm). Loops so low in the solar atmosphere (ef- 
fectively embedded in the chromosphere) cannot be visible in 
the Ly-g line; mo r eover, if compared with the observations of 
iPatsourakos et al. I (l2007h . our loops typicaUy have smaller pres- 
sures with respect to their estimates. 

Indeed, the shape of the cool loops we have found could give 
information on the orientation of the magnetic field in the tran- 
sition region, and in particular on its inclination. Since, as al- 
ready pointed out, these loops are very shallow, the magnetic 
flux tube should emerge at transition region temperatures with a 
quite small angle with respect to the Sun's surface and observa- 
tions should reveal a predominance of horizontal magnetic field 
direction at their height levels, assuming that these loops are al- 
most everywhere. 

We have also obtained quasi-static loops with maximum 
temperature in the range 2 x 10^ - 10'' K, using a realistic ra- 
diative losses function. These loops are smaller with respect to 
coronal loops but have diffe rent characterist ics compared to the 
static cool loops proposed bv lAntiochos & Noci (1986) and oth- 
ers. These loops in principle could be observed with current tele- 
scopes, but in order to resolve them in all their temperature ex- 
tension, we would need multi-temperature observations, i.e. dif- 
ferent UV lines formed at temperatures between 10^ - lO*" K 
with resolution of at least 1". Loops 17 and 18 have the maxi- 
mum temperature ar ound 2.5 x 10^ K that is the upper limit of 
the region in which iFeldmanI (Il983h located the "UFS (unre- 
solved fine structures)". We find that these intermediate temper- 
ature loops follow the scaling laws for coronal loops contrary to 
results of previous works based on the observational data (e.g.. 
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lBrowniri996l) . The loops obtained have indeed small pressures 
that make their parameters obey the RTV scaling laws, but these 
pressures are 1-2 orders of ma gnitudes low er than the ones esti- 
mated from the observations fBrownlll996h . Observed interme- 
diate temperature loops are usually associated with flows while 
we are studying possible solutions to the hydrodynamic equa- 
tions in quasi-static conditions without considering their role. 
The small values of the pressure we have found for these loops 
are the only possible in order to find solutions to the hydrody- 
namic equations in quasi-static conditions. Due to their small 
pressures, these loops obey to the scaling laws for coronal loops. 
The presence of substantial flows could influence their pressures. 

We are not able to reproduce the DEMs derived from obser- 
vations with only one set of parameters (a single lo op). A similar 
prob lem has already been reported in the paper of 'Susin o et al. I 
(12010) but for coronal loops. We find instead that a combination 
of cool and intermediate temperature loops, in particular pre- 
cisely due to their computed pressures, can give a DEM with 
a shape not too far from observed. This of course does not pre- 
clude the possibility that some additional physical processes, not 
included in our simulations (like that proposed by Judge 2008 
or bv iKlimchuk et al. I i2008 and De Pontieu et al. 1 12009. or dy- 
namic structures like cooling, heating loops and spicules), must 
be taking place in the transition region and contribute to explain 
the discrepancies with the observed DEM structure. Due to their 
dimensions, the emission of the quasi-static loops found can be 
seen as a diffuse component of the transition region. 

In our simulations, we have considered only the case of 
constant heating rate and we have already found quasi-static 
cool loops in conditions not allowed from strictly static and sta- 
tionary hydro-dynamical equations. More realistic assumptions 
could make obtaining stable, quasi-static cool loops even eas- 
ier In particular, in this work we have shown how the shape of 
the radiative losses function below 10^ K is very important for 
the existence of cool loops. Because the radiative losses around 
logT ~ 4.2 K are dominated by the H i Ly-a, we plan to explore 
in a follow-up study the effect on the structure and stability of 
cool loops of a more realistic treatment of hydrogen radiation 
losses in the lower TR. To do so, it will be necessary to improve 
the models by introducing a full treatment of optically thick ra- 
diative losses and/or partial ionization, and therefore solving the 
radiative transfer equations and rate equations of hydrogen (non- 
local thermodynamic equilibrium). Our first attempt to simulate 
an optically thick radiative losses function by removing the hy- 
drogen losses from a realistic radiative losses function seems to 
be promising. Without the hydrogen losses, the function Aoea-H 
becomes lower and we are able to obtain quasi-static cool loops 
with r < 10^ K. 

Taking into account finite cross sections for these loops could 
influence the details of the radiative transfer in the Ly-ff line (by 
reducing, for instance, the escape probability of Ly-o' photons). 
Predicting the eff'ect of loop area expansion above the chromo- 
sphere is not so easy, although it is difficult to envision substan- 
tial expansion factors in such small loops. 

Finally, the maximum temperature reachable by cool loops 
is limited only by the second maximum (or a change of slope) 
present in the A(T) since we have shown that the peak due to 
the Ly-a losses is not a limit. The different radiative losses func- 
tions used have this maximum (or the slope change) at different 
temperatures. So, we underline the importance of knowing accu- 
rately the A(r) also at temperature around 10^ K for the structure 
of the cool loops and of what we termed "intermediate loops", in 
analogy with the importance of the position and strength of the 
Ly-a peak for cool loops. 
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